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The design and operation of Beam- 
con III, the latest linearity measure- 
ment instrument using the beam 
addition method in the detector metrol- 
ogy program at the National Institute 
of Standards and Technology, is de- 
scribed. The primary improvements in 
this instrument are the reduction of 
stray radiation to extremely low levels 
by using three well-baffled chambers, a 
larger dynamic range, and an additional 
source entrance port. A polynomial re- 
sponse function is determined from the 



data obtained by this instrument using 
a least-squares method. The linearity of 
a silicon photodiode-amplifier detector 
system was determined to be within 
0.054 % (2<7 estimate) over nine 
decades of signal. 
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1. Introduction 

The underlying assumption in the use of most 
detection methods for the measurement of optical 
radiation is that the detector output signal is di- 
rectly proportional to the input radiation flux. This 
proportionality is defined as linearity and con- 
versely the departure from proportionality is de- 
fined as nonlinearity. Optical radiation detectors 
are typically converters of optical power (or photon 
flux) to measurable electrical parameters, such as 
current, voltage, or pulse frequency. Therefore, in 
radiometric applications, the linearity of a detec- 
tor's signal is not only a function of the efficiency 
of the physical radiation detector (transducer), but 
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also of the peripheral electronics. A well known 
example of this is a photomultiplier/photon-count- 
ing system where significant nonlinearity is induced 
at high count rates by the dead time of the pulse 
counter. All detector systems, consisting of detec- 
tors and peripheral electronics, are characterized 
by some degree of nonlinearity. The exact positions 
of these nonlinear regions are functions of both the 
detector and the accompanying electronics and 
should be determined experimentally for measure- 
ments of the highest accuracy. There are essentially 
two practical choices when confronted by nonlinear 
behavior of a detector system: first, to view it as a 
performance limit, a maximum deviation from lin- 
earity over a given range of input; second, to 
correct the measured signal given the known 
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nonlinearity of the detector. This latter solution 
presupposes that the detector's nonlinearity is well 
characterized and should be determined with the 
measurement electronics that are to be routinely 
used with the detector. 

Concern about detector linearity and the devel- 
opment of instruments for its measurement has a 
long history at the National Institute of Standards 
and Technology, formerly the National Bureau of 
Standards (NBS). Sanders [1] reported in 1972 the 
results of photocell linearity measurements using a 
multi-aperture design. Also in 1972, Mielenz and 
Eckerle [2] utilized a double-aperture method to 
characterize the NBS Reference Spectrophotome- 
ter. In 1984 Saunders and Shumaker [3] reported 
measurements using an automated radiometric lin- 
earity tester. This instrument was called a beam- 
conjoiner. 

Based on this initial work, we now report on our 
latest generation of beamconjoiner technology, 
Beamcon III. This instrument was specifically de- 
signed to characterize radiometric detectors, par- 
ticularly at low flux levels approaching the noise 
equivalent power (NEP) of the detector. To facili- 
tate low level measurements, this instrument has 
three well-baffled chambers to reduce the stray ra- 
diation to extremely low levels. Other improve- 
ments are additional filters to give a larger dynamic 
range and a second source entrance port for a laser 
source or for source mixing. The data analysis al- 
gorithm of Beamcon III calculates the responsivity 
of the detector system assuming a polynomial rela- 
tionship between the input flux and the output sig- 
nal from the detector system. Beamcon III will 
serve as an automated instrument for routine de- 
tector linearity characterization and for determin- 
ing the linearity of NIST's monochromator-based 
systems (e.g., spectroradiometers and spectropho- 
tometers). 



2. Description of Beamcon III and its 
Operation 

A schematic diagram illustrating the technique 
of beamconjoining is shown in Fig. 1. The input 
beam from the source is split into two optical paths 
by the first beamsplitter BSl and variably attenu- 
ated by filters on either wheel Wl or W3. After 
reflection off mirrors Ml and M2, the two beams 
are combined by beamsplitting BS2 and attenuated 
again by a third filter on wheel W2 before falling 
on the detector. All filters are metallized neutral 
density filters with quartz substrates, and the filter 
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Fig. 1. Schematic optical diagram of Beamcon III showing the 
primary optical elements. BS: beam splitter; W: filter wheel; 
and M: mirror. 

wheels are tilted to preclude interreflection be- 
tween filters. In addition the beamsplitters are 
wedged to prevent etaloning. 

A complete schematic diagram of Beamcon III, 
including all ports, optics, and chambers, is shown 
in Fig. 2. Beamcon III provides four basic advan- 
tages over the previous technology: 1) each cham- 
ber is optically isolated from the next chamber 
except for a 3 inch diameter throughput port, 
thereby reducing scattered radiation; 2) there are 
two source input ports for both laser and broad- 
band sources either singly or concurrently for 



Chamber 



d Detector 




Fig. 2. System layout of Beamcon III. Chamber 1 collimates the 
beam from the source. Chamber 2 splits the beam into two 
paths and recombines it. Chamber 3 has a spatial filter SF to 
eliminate stray light and focuses the beam onto the detector. 
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source mixing; 3) a steel base plate allows for mag- 
netic mounting of additional optical components; 
4) an external translation stage on the broadband 
input port can vary the numerical aperture and col- 
limate the source beam. A microcomputer controls 
the positions of the filter wheels and records the 
signal from the detector. 

An incoherent light source (Source, Fig. 2), ei- 
ther a quartz halogen or arc lamp, is placed on a 
translation stage in front of the input port to 
Chamber 1. The input power of the source beam is 
varied either by changing the aperture of the iris 
diaphragm (Iris 1) or by translating the Source, or 
both. The input beam in Chamber 1 is folded by a 
flat mirror (Mf) and collimated by a concave spher- 
ical mirror (Mc), which projects the beam through 
the input port of Chamber 2. For narrowband op- 
eration with an incoherent light source, an interfer- 
ence filter can be placed at the input port of 
Chamber 2. The first beamsplitter (BSl) yields two 
separate beams which are independently attenu- 
ated by two filter wheels (Wl) or (W3). These fil- 
ter wheels contain four neutral density filters with 
differing optical densities. The two beams are 
folded by mirrors (Ml) or (M2) and recombined 
into a single beam by beamsplitter (BS2). This 
beam is further attenuated by passing through a 
third filter wheel (W2) containing five different 
neutral density filters. The beam then passes into 
Chamber 3, which contains a set of concave mirrors 
(M3 and M4). A spatial filter (SF) at the focal 
point between the two mirrors filters out scattered 
and diffracted light, resulting in a uniform output 
field at the detector. A concave mirror (M5) fo- 
cuses the beam onto the detector. Translational 
shifting of M5 controls the image plane of the in- 
strument and varies the irradiance at the detector. 
The detector is generally overfilled to prevent arti- 
facts due to different cut-off edges in the optical 
path (optical stop variation) and to the local 
nonuniformity of the detector. In addition environ- 
mental parameters such as humidity and tempera- 
ture are constantly monitored at the detector. A 
monochromator can also be attached at the exit 
port to select a wavelength for narrowband opera- 
tion. 

Beamcon III has a second input port into Cham- 
ber 2, which can be used singly or in concert with 
the beam from Chamber 1. Generally the input 
port of Chamber 2 is used for sources such as He- 
Ne or Ar* lasers. A beam expander BE is used to 
minimize scattering losses and filter transmittance 
nonuniformity. The laser can also be used to align 
both optical paths by removing the reflecting mir- 



ror (Mr) and inserting the mirror (Ma) to guide the 
laser beam through the center of the lamp filament 
and of every optical element. The mixed source 
method has not been attempted as yet, but is 
planned in future studies, as it has advantages for 
studying differential spectral detector nonlineari- 
ties [4]. 

Stepping motors control the filter wheel posi- 
tions, each with 400 steps/revolution. A tolerable 
error is one step, therefore the error in filter posi- 
tion is 0.9°. The space between two filters on each 
wheel acts as a shutter position for a dark signal 
measurement. Thus, the total number of filter com- 
binations is 150, and a full 150 measurements is 
called a pass. Of these, 30 are dark signal measure- 
ments: 5 when both Wl and W3 are shuttered and 
25 when W2 is shuttered. The remaining 120 filter 
combinations, yielding 120 response signal mea- 
surements, are randomized in each pass to prevent 
hysteresis. The dark signal/filter combinations are 
randomized and inserted after every five measure- 
ments. Each dark signal is used to correct the suc- 
ceeding five response signal measurements. 
Generally, five passes are taken and averaged to 
reduce the random error and to collect information 
about possible source, filter, and detector drift. 
The maximum throughput is 7 % at a wavelength 
of 632.8 nm, and the total attenuation factor is 
about 3000 against a light source. 

The output beam is partially vertically polarized 
for an unpolarized source since each tilted surface 
preferentially reflects the vertical component of 
the incident beam. The ratio of the polarization in 
the vertical direction to that in the horizontal di- 
rection is 1.17. Thus, it is critical to keep the polar- 
ization direction of a source (e.g., laser) constant in 
order to eliminate this polarization preference. 



3. Polynomial Response Function Fitting 

The principle of operation of Beamcon III is the 
beam addition method, in which different filter 
wheel combinations result in different fluxes at the 
detector. This method assumes that the flux from 
the source remains constant during all measure- 
ments and that the output flux at the detector is 
the sum of the two fluxes along the two different 
paths. 

The final flux <f> at the detector is the sum of the 
fluxes <}>i and <f>2 from the two paths. Indexing the 
fluxes by the filter combination. 



<i>{ijJi) = 4>i(iJ<^) + <kiJJ<^\ 



(1) 
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where i equals to 4 designates the filter on wheel 
Wl,y equals to 4 the filter on W3, and k equals 
to 5 the filter on W2. In each case, an index value 
of corresponds to an opaque, shuttered position. 
For a single path, there are 30 possible filter com- 
binations. Of these, no flux is transmitted for 10 
combinations. Thus, each path has 20 filter combi- 
nations that result in a response signal. 

Since the detector is assumed to be nearly linear, 
it is appropriate and convenient to assume an n -th 
order polynomial response function for the detec- 
tor [3], so that 

(t>(ijJc) = ro + ris(ijJc)+r^^(iJJc) + ...+r„s"(ijJc), 

(2) 

y/here s (i j Jc) is the measured response signal, cor- 
rected by the dark signal and indexed by the filter 

combination, and ro, n, rz r„ are the response 

function coefficients. The coefficient n is set equal 
to one, both because the detector is assumed to be 
nearly linear and to normalize the flux to the same 
units as the signal. 

Combining Eqs. (1) and (2), with ri = l, yields 

+ ...+r„s"(ijJc). (3) 

There are 120 expressions of the form of Eq. (3) 
from the 120 distinct response signals measured in 
each pass, with 40+/i unknown variables from the 
20 filter combinations along each path and the n 
response function coefficients. The fluxes of shut- 
tered filter combinations are set equal to zero. 
These 120 expressions are solved using the linear 
least-squares technique to determine the values of 
the unknowns. While the values of the polynomial 
response function coefficients are the primary goal 
of the measurements, the values of the fluxes are 
useful for monitoring the stability of the filters. 

Solving for both the fluxes along the paths and 
the response function coefficients provides two 
ways in which to determine the total flux at the 
detector, using Eqs. (1) and (2), and therefore 
serves to verify the results. A measurement run 
yields 120 residuals between the total fluxes calcu- 
lated in these two ways. If there is a slope in a 
time-ordered plot of these residuals, then either 
dark signal or source intensity variations likely ex- 
ist. Also, the standard deviation of these residuals 
indicates the noise level present during the mea- 
surement, either from the source, filters, or detec- 
tor. For most detectors, a second-order polynomial 



accurately fits the data. However, if the detector is 
highly nonlinear, the power of the polynomial is 
increased until the residuals are normally dis- 
tributed about zero. If the response function coeffi- 
cients remain unreasonable, then either a 
polynomial is not an appropriate fitting function or 
the detector and its associated electronics are not 
functioning properly. 

The algorithm for determining the 40+/1 un- 
known variables from a complete measurement was 
extensively tested using synthetic data sets with 
known response function coefficients for second- 
order polynomials. It was also tested using syn- 
thetic data sets corresponding to two cases of 
nonlinear behavior, supraresponsivity and satura- 
tion at high flux levels [5]. The algorithm was suc- 
cessful in fitting these two behaviors with a third- 
and a sixth-order polynomial, respectively. 



4. Linearity Measurement of a Silicon 
Photodiode 

Beamcon III was used to determine the linearity 
of a silicon photodiode-amplifier detector system 
designed and constructed at NIST for multi-decade 
performance and described previously [6]. The lin- 
earity was determined for amplifier gains from 10^ 
to 10"". The broadband source was a 1000 W tung- 
sten-halogen lamp operated at a constant current 
of 7.6 A. In order to cover the entire signal range 
of the detector system without changing the spec- 
tral distribution of the source, the flux from the 
lamp was varied either by changing the numerical 
aperture of Iris 1 in Fig. 2 or by placing neutral 
density filters in Chamber 1, or both. 

The relative responsivity of a detector system in- 
dicates deviations from linearity. Once the re- 
sponse function coefficients have been determined, 
the total normalized flux for a given signal is given 
by the response function [the right-hand-side of 
Eq. (3)]. If the detector system is linear, the total 
normalized flux is simply the signal. Therefore, the 
relative responsivity at a given signal is the re- 
sponse function divided by the signal. The results 
from one measurement of the linearity of the de- 
tector system at a gain of 10* is shown is Fig. 3, 
where the relative responsivity is plotted as a 
function of signal. The response function was a 
second-order polynomial with coefficients 
ro= -5.53x10"' and r2= -1.56x 10-^ Over the 
three decades of signal measured at this gain, the 
detector system is linear to within 0.027 % (2o- 
estimate). 
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Fig. 3. Relative responsivity as a function of signal for a silicon 
diode-amplifier detector system at a gain of lO*. Tlie error bar 
indicates tlie 2<t estimated uncertainty of tlie relative responsiv- 
ity. 

The maximum signal for the detector system, de- 
termined by the source and throughput of Beam- 
con III, was about 10"' A. At a bandwidth of 0.67 
Hz, the NEP of the detector system resulted in a 
minimum measurable signal of about 10"'^ A. In 
order to characterize the detector system between 
these limits, linearity measurements were per- 
formed at all seven amplifier gain ranges. A sec- 
ond-order polynomial response function was 
calculated for each range, and the relative respon- 
sivity over nine decades of signal is shown in Fig. 4. 
The detector system is linear to within 0.054 % (2o- 
estimate) over the entire range, and to within 
0.209 % (2o- estimate) for currents less than 10" " 
A. Thus, Beamcon III was able to determine the 
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Fig. 4. Relative responsivity as a function of signal for a silicon 
diode-amplifier detector system over its entire signal range of 
nine decades. The error bars indicate the estimated 2a- uncer- 
tainty of the relative responsivity for signals less than 10" A (on 
the left) and over the entire range of signal (on the right). 



linearity of this detector system for signals ap- 
proaching its NEP. 

5. Conclusions 

The design and operation of our latest genera- 
tion of beamconjoining technology, Beamcon III, 
was detailed. Using this instrument, the linearity of 
a silicon photodiode-amplifier detector system was 
determined over nine decades of signal, approach- 
ing its NEP at the smallest signals. Beamcon III 
will permit both routine calibration of detector sys- 
tems and research into detector-amplifier effects 
on nonlinearity. Future studies on the linearity of 
solid state detector-amplifier systems, photomulti- 
plier photon-counting systems, and multichannel 
detectors (i.e., diode arrays and charge coupled 
devices) are planned, as well as on the capabilities 
of Beamcon III for source mixing. 
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